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We investigate a helicity non-fiip double transverse spin asymmetry in vector boson pro- 
duction in hadron-hadron scattering, which was first considered by Ralston and Soper at the 
\ tree level. It does not involve transversity functions and in principle also arises in W-boson 

. production for which we present the expressions. The asymmetry requires observing the trans- 

' verse momentum of the vector boson, but it is not suppressed by explicit inverse powers of 

, a large energy scale. However, as we will show, inclusion of Sudakov factors causes suppres- 

^ • sion of the asymmetry, which increases with energy. Moreover, the asymmetry is shown to be 

' approximately proportional to Xigi{xi) X2gi{x2), which gives rise to additional suppression at 

, small values of the light cone momentum fractions. This implies that it is negligible for Z or 

W production and is mainly of interest for 7* at low energies. We also compare the asymmetry 
with other types of double transverse spin asymmetries and discuss how to disentangle them. 
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■ I. INTRODUCTION 

o ■ 

I Double transverse spin asymmetries in high energy hadron-hadron collisions have attracted much theoretical 
attention (starting with the early investigations but no experimental studies have been performed so 

far. At the BNL Relativistic Heavy Ion Collider (RHIC) transversely polarized protons will be collided for 
^pH, the first time. Therefore it is important to make an analysis of the double transverse spin asymmetries. The 
O 1' proposed DESY HERA- TV experiment also prompts such a study (see e.g. |^). For the transversity double 
O I spin asymmetry in the Drell-Yan process such studies have been performed in considerable detail [Q. In this 
1^ ■ paper we will mainly investigate another type of double transverse spin asymmetry, one which does not involve 
transversity functions. 

In general, helicity non-flip quark and gluon states in transversely polarized hadron-hadron scattering will 
^ . lead to power suppression [0{1/Q'^), where is the vector boson virtuality]. In the present study we will 

■ exploit the fact that if the transverse momentum of the produced vector boson is observed, this no longer holds 
true. By observing the transverse momentum certain azimuthal asymmetries can occur in the cross section 
without explicit power suppression. The possibility of such an unsuppressed helicity non-flip double transverse 
spin asymmetry has been noted before in the literature and a tree level expression for the case of a virtual 
photon has been given |jl|,|| . Here we investigate speciflcally the case of weak vector boson production and the 
effect of inclusion of Sudakov form factors. 

The helicity non-flip nature will allow for double transverse spin asymmetries even in W production (for which 
the helicity flip contribution is absent Q). Unfortunately, such transverse momentum dependent azimuthal 
asymmetries will turn out to suffer from suppression due to Sudakov factors, which increases with energy, as 
was forseen in Ref . 0| . We will explore this issue quantitatively in detail and we will show that for the azimuthal 
double transverse spin asymmetry of interest, the inclusion of Sudakov factors causes suppression by at least an 
order of magnitude compared to the tree level result and effectively produces a power behavior of l/Q", with 
a ~ 0.6. Moreover, the asymmetry will be shown to be approximately proportional to a;igi(a;i) 2:232 (2^2), which 
gives rise to additional suppression at small values of the light cone momentum fractions. The conclusion will 
be that this asymmetry is of interest mainly at lower energies, i.e. for 7* production. This also leaves the option 
of studying possible contributions to double transverse spin asymmetries in W production from physics beyond 
the standard model. 

The outline of this paper is as follows. In Sec. II we will repeat the essentials of the transversity double spin 
asymmetry in the Drell-Yan process in order to contrast it to the helicity non-fiip asymmetry (Sec. HI). We 
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study the latter asymmetry in the neutral (Sec. IIIA) and in the charged (Sec. IIIB) vector boson case. In order 
to obtain estimates we will assume Gaussian transverse momentum dependence of the quarks (discussed in Sec. 
IIIC). We will then include Sudakov form factors in the asymmetry (Sec. HID) and estimate its quantitative 
effects (Sec. HIE). In Sec. IIIF we will comment on the possible study of physics beyond the standard model 
via double transverse spin asymmetries in W production. 



II. TRANSVERSITY DOUBLE SPIN ASYMMETRY 



The main characteristic of the transversity double transverse spin asymmetry of vector boson production 
in hadron-hadron collisions is that the gluon distribution does not contribute. Hence, at leading order in an 
expansion in inverse powers of the hard scale(s) only the quark transversity distribution function (denoted 
as hi, Sq or Atq) contributes. This leads to the well-known expression for the double transverse spin asymmetry 
in the Drell-Yan process: 



a{p^ ££' X)-a{p^ ^ ££' X) _ sin^ 61 cos Eg.a Sg ^^^2) 
a{pf p^ ^ ££' X)+a{p^ pi- ££' X) ^ 1 + cos2 6I T^aa^lfili 



(1) 



We would like to note that the sum is over flavors including anti- flavors, otherwise one should add a term in 
both numerator and denominator for the exchange (2:1 ^2:2), since one can use that hi — h^. The above 
asymmetry comes from the following azimuthal dependence in the cross section 



dajp^ p^ ££' X) _ 



Y.el{y{l~y)\S,^\\S^T\ cos(4^ + 4j /i? (xiK(x2) + . . . }. (2) 
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The above is expressed in the so-called CoUins-Soper frame [|j (for details see e.g. [ p^pT| ): 

i=q/Q, (3) 

Xl ~ X2 ~ ... 
h = qT/QT^iq-XiPl-X2P2)/QT, (5) 

where q is the vector boson momentum, Pi are the hadron momenta and Pi = Pi^q/ {2xi). The azimuthal angles 
lie inside the plane orthogonal to t and z. In particular, (f) gives the orientation of = (gf"" - i'^i" + z^z") l^, 
the perpendicular part of the lepton momentum I; (j)g, is the angle between l± and Sjj^. In the cross sections we 
also encounter a dependence on y = l~ /q~ , which in the lepton center of mass frame equals y = {1 + cos0)/2, 
where 9 is the angle of z with respect to the momentum of the outgoing lepton /. 

The perturbative corrections to the double transverse spin asymmetry, Eq. (|^), have been calculated in fl^ 
and using the assumption that at low energies the transversity distribution function hi equals the helicity 
distribution function gi or by saturating the Soffer bound, it has been shown in Ref. that Att is expected 
to be of the order of a percent at RHIC energies. 

In Ref. |q| it was discussed that in the process p^ p^ W X the transversity distribution cannot contribute 



and this is a general feature of chiral-odd functions and charged current exchange processes (see also |13 ). 
This means that only the suppressed contributions from the twist-3 distribution function gx and its gluon 
analogue Axg contribute (these are chiral-even functions; they mix under evolution). Of course there are 
contributions of the transversity functions via quark mass terms or via production of other particles that can 
compensate for the helicity flip, but these are all of higher order in the strong and/or weak coupling constants 
(e.g. one can think oi p^ p^ Z X W~ X , but this is negligible at RHIC energies). Since neither 

quarks nor gluons contribute without suppression to the asymmetry Att in the process p^ p^ — > W X, it might 
make this asymmetry a good place to look for contributions from physics beyond the standard model. For 
instance, scalar or tensor couplings of the quarks to the W could in principle produce an asymmetry. We will 
return to this issue at the end of this article. First one has to investigate and estimate another standard model 
mechanism, namely, there is the possibility that the quarks (and gluons) are not exactly collinear to the initial 
proton, leading to a helicity non-flip asymmetry without explicit suppression factors of 1/Q. In other words, if 
one measures the cross section differential in the transverse momentum of the vector boson, either in its angle 
compared to the other particles or in its magnitude, the helicity non-flip double transverse spin asymmetry can 
receive contributions at leading order, even for W production. If one averages over this transverse momentum, 
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then the asymmetry will vanish, but an (inadvertent) incomplete averaging, for instance due to imposed cuts, 
might still have observable consequences, cf. for instance Ref. Even though we will show that for W 

production this will not be a problem since the asymmetry turns out to be negligible, for 7* at lower energies 
this is important to take into consideration. 



III. HELICITY NON-FLIP DOUBLE TRANSVERSE SPIN ASYMMETRY 



If one can measure the cross section differential in the transverse momentum of the vector boson, either in 
its angle compared to the other particles or in its magnitude, then there is the possibility to have a double 
transverse spin asymmetry at leading order, in principle even for W production. To illustrate this we will 
make use of the formalism pioneered by Ralston and Soper which will be applicable in the region where 
the observed transverse momentum is small compared to the hard scale(s). Since this is a tree level formalism, 
we will later on include the effects of resummation of soft gluons by combining it with the approach of Ref. 
Ip^ . We focus on the Drell-Yan process, first on neutral vector boson production and later on charged vector 
boson production. However, the expressions also apply to p^p^ — *■ 2 jets (same formulas with minor trivial 
replacements) and a similar analysis can be applied to ^ tt (A, . . .)X, where the observed hadron is part 
of a high-pT jet. But in those cases there are background processes, which should be considered also and this 
requires more detailed study. 



A. Neutral vector boson production 

In Eq. (|^) we have given the contribution to the cross section that depends on the sum of the two transverse 
spin angles with respect to the lepton pair production plane, i.e. cos(0g^ + (f>g^)- This means that if one 
integrates over the lepton pair orientation, then this azimuthal dependence will average to zero. At order l/Q^ 
both quarks and gluons can contribute to Att via a term in the cross section which does not depend on the 
lepton scattering plane 

Att oc cos{(j)s^ - (j)s.^ ) gr 9t, (6) 

but this is expected to be negligible at — M|. Moreover, it is not at all clear that such a factorized 
description of the asymmetry holds at this level of next-to- next-to-leading twist, since it is well known that for 
the unpolarized case this order [©(l/Q^)] does not factorize. 

On the other hand, if one were to observe the transverse momentum of the lepton pair compared to the 
protons, there will be a double transverse spin asymmetry as a function of this transverse momentum qt, which 
appears at the leading order in 1/Q. It will involve one more angle (0^), but even if one would integrate over 
this angle (keeping only the magnitude of qt) and over the lepton pair orientation (0^), then there will remain 
an azimuthal dependence in the cross section that depends on the orientations of the two transverse polarization 
vectors only. 

To make this explicit we will look at Eq. (Al) of Ref. [0, which gives the cross section for the polarized 
Drell-Yan process p"^ p^ ~* 7 {Z) X in the formalism of Ralston and Soper using transverse momentum 
dependent distribution functions. From the expressions for the production of the Z boson it is easy to obtain 
the expressions for the production of the W boson. We will not repeat all the details of the calculation of the 
cross section expressions, rather we will focus on the expression as given in the Appendix in Ref. [ pT| . For 
the sake of argument it is unimportant to include contributions from the (formally difhcult) T-odd distribution 
functions, hence we neglect them, but they can be easily included. This leaves 



d.{h,h, ^ IIX) = ^ ^ J j^a^y^ \s,t\ \S,t\ sH4 - A) Bin(/. - 4J ^ 



d^dxidx2<P Qt 

-K^iy) \S,t\ I^2tI cos(2^1 - 4^ - 4j T 

where 



h prp h krp _ 



(7) 
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-y + y^) [el + 2g'yeag^xi + c[c1x2] 



Ki{y) 



which contain the combinations of the couphngs 
The Z-boson propagator factors are given by 



l-2y 



X2 



1 Q2(Q2 _ Ml) 



1 

Xi, 



(8) 

(9) 
(10) 

(11) 
(12) 



and gv and gA are the vector and axial- vector couphngs to the Z boson. We have summed over the polarization 
of the outgoing leptons. Furthermore, we use the convolution notation (Ralston and Soper |^ use /[...]) 



J^[ff] = J <fpTd?kT5\pT + kT-qT)nXl,pl)r{x2,kl), 



(13) 



where a is the flavor index. 

The function giT is the function h^^ of Ralston and Sopcr and has as interpretation the distribution 
of longitudinally polarized quarks (7^75 projection) inside a transversely polarized hadron. It enters into the 
calculation compared to the unpolarized distribution function as follows: 



$(xi,pt) 



Ml 
2P+ 



h{xi,PT) 



2 ^ 



Ml 



{Pt ■ Sit) 

Ml 



9it{xi,Pt) 



ri75 

Ml 



(14) 



The details of the momenta are: the momenta of the quarks which annihilate into the photon with momentum 
g, are predominantly along the direction of the parent hadrons. One hadron momentum (Pi) is chosen to be 
along the lightlike direction given by the vector n+ (apart from mass corrections). The second hadron with 
momentum P2 is predominantly in the n_ direction which satisfies n+ -ri- = 1, such that Pi ■ P2 — 0{Q^). We 
decompose the momenta in +, — and transverse components, defined through = p ■ nip, where we note that 
[cf. Eqs. i)-(|)] 



1 

1 
72 



Q 



Q 



(15) 
(16) 



The four-momentum conservation delta function at the vector boson vertex is written as (neglecting 
contributions) 



S^q -k-p) = 5{q+ - p+) 6iq- - k') S^{pj. + - Qt) 



(17) 



and allows for integration over p~ and fc+. However, the transverse momentum integrations cannot be separated, 
unless one integrates over the transverse momentum of the vector boson or -equivalently- of the lepton pair. 

The Drell-Yan cross section is obtained by contracting the lepton tensor with the hadron tensor. At the tree 
level we find, for the hadron tensor. 



W^'^ = - J d^pr<fkrS^{PT + kT^qT) Ti-mxi,Pr)Vr<i>{x2,kr)V2%+^_ + i^'^^^^jy (18) 

The vertices Vf can be either the photon vertex V = 67^^ or the Z-boson vertex = gvl^ + gAlbl^ ■ 

The above given azimuthal dependence in the cross section, Eq. (^, means that if one observes the transverse 
momentum of the 7 or Z boson, one can consider the cross section differential in the magnitude of the transverse 
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momentum only and integrate over the orientations of the leptons and of qx itself. This results in the following 
double transverse spin asymmetry: 



da [p^ ^ ^ jZ) X] -da[p^p^ ^^{Z)X] _ Y.a,-a Kfiy) T [p^-k^ girg^r] 
^^^"^^^ da[p^p^ ^^{Z)X]+daWpi^^{Z)X] 2M,M, ^^-^ K^{y) T [hJ,]' ^ ' 

where Q\ = —q^ = q"^ <^ Q"^. This can be seen from the following considerations. The angular dependence 
sm{(j)f^ — (j)gj s'm{(j>j^ — (j>g^) can be rewritten after integration over the angle 0^ as cos{(j>g^—(j)g^)/2 = cos((/)|j)/2. 
Since this does not depend on the orientation of £ itself, one can integrate over it also. The angular dependence 
cos(20^ — — averages out. 

If on the other hand one only observes the angle of the transverse momentum and averages over its mag- 
nitude, one can also obtain a nonvanishing asymmetry (which can still be differentiated from the transversity 
asymmetry). The whole point is to prevent the averaging: J d'^q^^ !F [prp-kj^ giTg^rp] = 0. As said before, an 
incomplete averaging due to imposed experimental cuts might also result in a nonvanishing asymmetry, cf. for 
instance Ref. p^ ]. 

Before we continue we would like to point out that unlike for hi there is a leading twist gluon analogue of the 
function giT- The function arises with ie^^ in the correlation function cx {PS\F^°'F^^\PS), which means 
that it is a At/ type of function with transverse momentum dependence. But since the transition g g is 
only possible via quarks, we implicitly include the gluon in the sum over flavors. 



B. Charged vector boson production 

In order to arrive at the expressions for the cross sections of the charged current process (cf. Ref. one 
can take = and replace 

X2 X2 ^ (o.:„2/, ) 77t2 'nTyTT'FmT ' (20) 



V 6 sin^ Ow J iQ' - + r^M^, 

in the above given coupling if". In addition, one replaces ci = ±C3 — 1, depending on the chirality of the quark 
or lepton, since ci — {gj^ + g'j^)/'2 and C3 = {gj^ — gjf)/^. Hence, for a left-handed quark one finds c° = = 1 
and for a right-handed quark one finds c° = — — 1; similarly for the leptons. We also note that a left- or 
right-handed quark or lepton has helicity A^/g = Tl- This results in 



— 4 ^IV P X < ^'^^ equal quark and lepton chiralities, 

^'^ ' I (1 _ y)2 opposite quark and lepton chiralities, ^ ' 

where a, b are the incoming quark and antiquark flavor indices, respectively, and Vab stands for the appropriate 
Cabibbo-Kobayashi-Maskawa (CKM) matrix element. We illustrate the above by assuming that only u and d 
quark distribution functions contribute. This leaves two elementary subprocesses: ud — > — s- e^iy {u and u 
have equal chiralities) and du — > W" — > e~P (d and e~ have equal chiralities) for which one finds the couplings 
K'^d- = K"'^ = ^x^\Vud\'^y'^ ■ For the cross section one has to take into account that in the sum over final state 
polarizations there is now only one state that contributes, but for the asymmetry this is not relevant. 

In case oip^ p^ X and subsequent leptonic decay of the W, we encounter the problem that the produced 

neutrino will prevent a determination of the transverse momentum of the W boson. Hence, in the case of a 
produced neutrino one cannot define a lepton scattering plane (one does not observe T), hence azimuthal angles 
cannot be defined compared to the W boson direction. This holds unless one can reconstruct the direction of 
the neutrino by the momentum imbalance [^|. 

Another possible way of observing the transverse momentum of the W boson is looking at the W decay into 
2 jets. The expressions for lepton pair production stay essentially the same for the 2 jets case after the obvious 
replacement of the coupling constants. By measuring the direction of the 2 jets, the transverse momentum of the 
W can be determined, but one problem is that it has ^* /Z 2 jets as a background. Separation of 'y* /Z and W 
might only be possible with a very high transverse momentum cut JTst , but then the given expressions are not 
applicable anymore. Another problem is that it also receives contributions from quark-quark scattering next to 
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the quark- antiquark scattering, but such contributions have rather large color factor suppression In any 

case, the contribution coming from the transversity functions to the 2 jets asymmetry can always be eliminated 
by averaging over the orientation of the 2 jets, as explained above. 

Here we will focus only on lepton pair production and assume that in case of W production the direction of 
the neutrino can be reconstructed to obtain the transverse momentum of the W. 



C. Gaussian transverse momentum dependence 



In order to obtain an estimate of the above asymmetry, we will consider a Gaussian transverse momentum 
dependence of the functions. Instead of using Gaussians, another way of obtaining an estimate of the asymmetry 
would be to use the spectator model for the function giT{x,p'^) 20 1 . But for simplicity we will assume a Gaussian 
transverse momentum dependence, e.g. 

9it{x,p'^t) = 9it{x) — exp {-R^pl) = gixix) Gip^r)- (22) 

We would like to relate the function giT{x) to a well-known function in order to be able to make some 
predictions in the end. This can be achieved by using p|, weighted functions 

(23) 

For a Gaussian transverse momentum dependence we find that 

g^iT giT{x)l{2A'PR^). In the Wandzura- 
Wilczek approximation the function gji^ is a well-known quantity: it equals (upon neglecting quark masses) 
xg^^ {x), where g^^ {x) = gi + 32^^ is the Wandzura-Wilczek part of the function gx (see also Ref. |2l| 
for a discussion on this topic). This can be shown by using the equations of motion. The function gT has 
been studied by SLAG and the Spin Muon GoUaboration (SMC) j2^ and the data are (still) consistent with 
9t = g^^ ■ Also, the data show that 52 is small compared to gi, therefore, up to a few percent one can take 
xg^^{x)Kixgi{x). Thus, we find giT{x) = g[l^2ApR^ ^ xgi{x)2]VpR^. 

Furthermore, will assume that the Gaussians are the same for /i and gir and for both protons, i.e., we take 
Ri = R2 ^ R and Mi = M2 = M. In this way we find for example 

T[pT-kj.giTgiT]~ ~[ 2~)^^'^\ 2~) ^^Siixi) X2gi{x2). (24) 

This results in the following tree level double transverse spin asymmetry at Qt = 

. nJ2nr, Kf(y) xigf(xi) X2gUx2) 
Att(Qt = 0) = M^R"" ^ — iv^y iyi\ ij jgiK 2) ^ 

Eara K<t{y) /f (Xi) n{x2) 

We will also encounter the Fourier transforms of these functions. The function / will denote the Fourier 
transform of /, and since we use the notation f{x) = J cPpj, f{x,prp), we see that f{x) = /(x, b = 0). Taking 
the Fourier transform of Eq. ( p^ yields 

g^j,{x, 5^) = giT{x) exp (-j^^ ■ (26) 



D. Beyond the range of intrinsic transverse momentum 

Monte Garlo studies including soft gluon resummation show that the largest contribution to the unpolar- 
ized cross section arises when the transverse momentum of the W is several GeV. This transverse momentum 
is too high to trust tree level expressions which involve only intrinsic transverse momenta. In order to go 
beyond this region, we will also include the Sudakov factor arising from resummed perturbative corrections to 
the transverse momentum distribution. 

Resummation of soft gluons into Sudakov form factors [p3 results in a replacement in Eqs. 



13) and (M) of 
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(27) 



where e ^^^"^ is the Sudakov form factor and — b^. This has been shown in Rcfs. p5| , p^ for the leading twist 
and is discussed for the present context in more detail in Ref. [p6|. The Sudakov form factor is found to be 



(28) 



One can expand the functions A and B in and the first few coefficients are known for unpolarized scattering 
prf and for longitudinally polarized scattering |2|]. The latter result is needed here since the function giT is a 
distribution of longitudinally polarized quarks; the asymmetry on the parton level is all- In order to obtain a 
first estimate of the effect of including the Sudakov factor we will take into account only the first term in the 
expansion of A: A^^^^ = Cf/tt- This leads to the expression 



S{b,Q)^ 



16 



33- 2n 



log 



log ( ^ ) log 



log(Q2/A2) 



(29) 



with bo = 2exp(— 7^;) « 1.123. We will take for the number of flavors nj = 5 and also Aqcd = 200 MeV. 
The replacement in Eq. (Oh leads to (suppressing the flavor index) 



f{xi,b) f{x2,b) 



— / dbbJo{bQT)e-^^''^f{x,,b)fix2,b). 

Jo 



(30) 



The functions also have a renormalization and factorization scale dependence, which we will choose to be equal 
fJ-R = fJ-F = Hence, we have, for instance, 



(fpj,f{x,pT\ij) = f{x,b = 0;/i), 



(31) 



where also the boundary of the integration gives a fi dependence. In Eq. ( |30| ) one usually takes f{xi, b; fi — bo/b) 

Of course, if one includes the effects of perturbative corrections, one should also include higher order cor- 
rections to the hard part not coming from soft gluons. But since the formalism, which means the factorized 
formula, is valid only for Q^, we include only the effects of soft gluons, which should allow us to extend 

the range of applicability from the region of intrinsic transverse momentum to the region of moderate Qrp val- 
ues. An equivalent way of saying this is that one can perform a collinear expansion of the hard scattering part 
H^'^ {xi,X2,PTikr]^,qj',Q) « H'^'^{xi,X2,Q), such that perturbative corrections to the hard part do not affect 
the transverse momentum structure of the tree level result and hence the transverse momentum weight in the 
asymmetry will be the same also beyond the tree level. Here we will consider H^'^{xi, X2, Q) to lowest order in 
as, therefore, only logarithmic corrections to the results presented below are expected. 

The numerator in Eq. (091) cannot be treated exactly like the denominator, so let us focus next on 



(32) 



T[pr-krff]^ J -0^ 6'" '^- e-^W J <fpT d^k^ pr ■ e-*(P-+'=-) f{x^,p'r)fi^2, fc|). 



As mentioned before, we assume that the distribution functions are Gaussians (as a function of transverse 
momentum), all of equal width: f{xi,p'^) = f{xi) G{Pt) ^nd f{x2,k'^) = f{x2) G{k^). One can then change 
variables to m = (pj. -I- kj,)/^/2 and v = {pj, — fcj,)/V2 and compute 



j d^pT d^k^ Pt ■ kj- e~'^<PT+kr) g{^pl) g{kl) = cxp 



2i?2 ) ' 



(33) 



which after application to Eq. (p2h yields (see also Eq. (M)) 
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r[PT-krff]=- j (0je*'^-^e-^^'V>i,fc)/(^2,6), (34) 

which can be compared with Eq. (|2^). Both equations fulfiU the property that the expression should vanish 
after integration. One infers that the numerator of Att{Qt) hence Att{Qt) itself oscillate, but in 
general / dQj, Att{Qt) 0; because of the Qj, dependence of the denominator. 
For the asymmetry, Eq. ([l9|), we then obtain 



Ea,-a Kt{y) X2g'i{x2) J,"° db MbQ t) exp [-5(6) - i^Vi?" 



2 Ea.-a Kt{y) m^l) mX2) C dbbUbQT) CXp [-5(6) - i b^ / R^] 



(35) 



where the approximation arises from taking giT{x) = gjj! 2M^R^ « a;(7i(a;) 2M^R^. 

In order to extend the above equation to the nonperturbative region of large values of 6, one usually introduces 
6* = 6/ yjl + 6^/6^3^^ and an additional term exp [— 5'Arp(6)], needed to describe the low qr region properly. In 
part, SNp{b) is introduced to take care of the smearing due to the intrinsic transverse momentum, therefore, 
taking into account the term exp (— i 6^/_R^) in addition will just produce a change in the coefficient of the 6^ 
term in SNp(b). To keep the unpolarized cross section unaffected, we will therefore introduce as nonperturbative 
term exp [—SNp{b) +5 b"^ / R^\ and study the following final expression for the asymmetry: 

Att{Qt) = - ^ - A{Qt), (36) 

2 Y.a,-a my) m^i) m^2) ^ ' 

where we define 

. ^ .,2 C dbb^ MbQr) exp [-5(6,) -SNp{b)] 
Ji[ljT)-M j-rf5^j^(;,Q^)exp[-5(6,)-5Arp(6)]- 

The denominator is then the conventional unpolarized expression. Also, we note that A{Qt) is dimensionless 
and, for simplicity, we will take M — 1 GeV. 

The above approach of including Sudakov factors in the tree level azimuthal asymmetry expressions can 
also be applied to expressions derived in for electron-positron annihilation and in ||l^ for lepton-hadron 

scattering. 



E. Estimating the asymmetry 

For the case of W production the asymmetry becomes 

A-^AQt) = - ^^^^^ Wa^? -^9l{^^) -29\{X2) 

If only the u and d quarks contribute, then also the CKM matrix elements drop out of the ratio. 

For the nonperturbative Sudakov factor we use the parameterization of Ladinsky-Yuan, Ref. p2{ , which was 
fitted to relevant Fermilab data, 

5Arp(6) = 5162 + 51536 ln(100xia;2)+.g262 In (^^^ , (39) 

with 51 = 0.1IGeV2,52 = 0.58GeV2,53 = -1.5GeV-i,Oo = 1.6 GeV and 6max = 0.5 GeV-^. We will take 
X1X2 = 10^2^ which is justified below. This leads to 5jvp(6) = 1.98 6^ at Q = 80 GeV and 5Arp(6) = 0.77 6^ at 
g = 10 GeV. 

The result for the asymmetry factor A{Qt) at Q = 80 GeV is given in Fig. |l|. 
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FIG. 1. The asymmetry factor A{Qt) a.t Q = 80 GeV. 

It is plotted up to Qt — Q/W, since beyond that Qt range the magnitude only slowly decreases and also 
the approximation Qt ^ Q^ is expected to become less valid. The asymmetry factor A{Qt) at Qt = is 
seen to be around 0.32 and at Qt values of a few GeV -relevant for the majority of produced W bosons- the 
asymmetry factor has a sign change and consequently a smaller magnitude. On top of that the asymmetry is 
proportional to \xigi{xi) X2gi{x2)\ < X2fi{x2)- Therefore, the total asymmetry as a function of Qt is 

expected to be below the percent level, if one assumes that on average xi = 0.4 and X2 = 0.07 for W production 
at RHIC ll^. Since we have implicitly included the gluons in the sum over flavors, the latter argument is not 
valid if Ag turns out to be extremely large at small x. Of course this will have even more serious implications 
for e.g. Aj^i^ in jet production, especially at low transverse momenta. We will just view this unlikely option as 
a proviso. 

In Ref. Q it is demonstrated that the transversity double spin asymmetry, which is a Qt integrated asymmetry, 
is expected to be at most on the order of a few percent, which matches the level of sensitivity of RHIC. The 
present asymmetry is still a function of Qt, requiring even more statistics. This will make the asymmetry 
^TT^'' (Qt) invisible at RHIC. Moreover, since it oscillates as a function of Qt, one expects that the asymmetry 
partly integrated over Qt, will not lead to any significant result either. 

A few remarks about the dependence of the result on the nonperturbative parameters. The asymmetry factor 
is seen to decrease with increasing Gaussian smearing width. Taking a higher value of &max and a lower value 
of xiX2 both increase this width. The above -optimistic- choices of &niax = 0.5GeV~^ and xiX2 = 10~^ are 
therefore expected to overestimate the asymmetry factor somewhat. 

But at lower energies -where larger light cone momentum fractions can be achieved- this asymmetry for 7* 
production is still worth investigating. 




1 2 3 4 5 

Qt [GeV] 

FIG. 2. The asymmetry factor A{Qt) at Q = 10 GeV. 
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In Fig. ^ we have given the function A{Qt) at the scale Q = 10 GeV and find that at low values of Qt it is 
around 1. Measuring A^j^(Qt) at larger values of xi and X2 might make this asymmetry observable. 

If one studies the values of A{Qt) at Qt = -where the asymmetry is largest- as a function of Q, one 
observes that inclusion of the Sudakov factor produces to good approximation a 1/Q" behavior, with a « 0.6. 
Even though this suppression is not very strong as a function of Q, one actually needs to compare the resulting 
asymmetry Eq. with the tree level asymmetry Eq. (25). This leads to a comparison of A{Qt = 0) and 
2M^i?^. In a tree level analysis = l/{p\) for a typical intrinsic transverse momentum squared value, i.e. 
R? Ki 2 — 1\ GeV~^, corresponding to the range of {p^) « (300 — 700 Me V)^. If we view i?^ = 1 as giving a 
lower bound for the tree level asymmetry factor, then A{Qt = 0,Q ~ 80) is still an order of magnitude smaller. 

We conclude that the Sudakov factors produce a strong suppression compared to the tree level result and 
the effect increases with energy. This will also have consequences for similar types of transverse momentum 
dependent azimuthal asymmetries appearing in for instance e~^e^ annihilation at the Z mass scale [ ^Tp^ , where 
the same strong suppression due to Sudakov factors is expected. 



F. Physics beyond the standard model 

It is now clear that the standard model (SM) mechanisms seem to produce negligible double transverse spin 
asymmetries in W production (A^rp). In summary the reasons are the following: the transversity distribution 
hi does not contribute At next-to-next-to-leading twist [0{MiM2/Q^)] the twist-3 distribution function gr 
(which is a helicity non-flip quark distribution) can contribute and its gluon analogue Axg as well. At = 
these contributions will be negligible. Furthermore, one can also neglect contributions which are of higher order 
in the strong and/or weak coupling constants. In the case of perturbative QCD corrections double helicity 
flip will be accompanied by quark mass terms and therefore will also be suppressed by at least two factors of 
1/Q- In the case of weak corrections one also expects negligible contributions, at least at RHIC energies, say 
around = Myy. Any leftover asymmetry from incomplete averaging of the transverse momentum dependent 
asymmetry A^t{Qt) we have estimated to be negligible as well. So at RHIC energies A^j. is expected to be 
negligible within the SM. 

It is now fair to address the question: if a significant asymmetry would nevertheless be found in the polarized 
proton-proton collisions at RHIC, can one really conclude something about physics beyond the SM? For instance, 
there could be scalar or tensor couplings of the W to quarks that can generate an asymmetry A^j,. If the scale 
of such new physics is A ^ then one might need to compare effects of order M\M2lQ'^ with / . 

For instance, for = and A = 1 TeV the latter is a factor of 40 larger (although still quite small). 

But the problem of comparing to higher twist contributions disappears altogether if the new couplings violate 
symmetries. 

This means that on top of the fact that the various SM mechanisms produce negligible asymmetries, one can 
also exploit the dependence on the orientation of the transverse spins compared to the lepton production plane 
to cancel out specific contributions exactly. We have seen that the transversity asymmetry ATJ^ appears with 
an angular dependence cos((/)g^ -f <t)%^): whereas Att{Qt) and the l/Q^ suppressed contribution from qt^t i'^ 
Eq. (^) both appear with cos(0g^ — The latter does not depend on the lepton scattering plane, because 

the asymmetries are not double transverse spin asymmetries at the parton level. On the other hand, symmetry 
violation asymmetries can produce other angular dependences than any SM mechanism. 

There might be T-odd asymmetries, for example the one of Ref. p^ , 

A^j, (X sin(4^ + 4j, (40) 

which would arise due to CP-violating vector couplings of the quarks to the W, which is assumed not to be V—A 
anymore, but some complex linear combination of V and A. It can clearly be distinguished from possible initial 
state interaction effects, which are P-even and only lead to asymmetries independent of the lepton scattering 
plane. To be more specific, if one assumes T-odd (P-even) distribution functions to be nonzero, then there will 
also be contributions proportional to [cf. Eq. (Al) of JlT| ] 

cos(/s^ - (fk )^[PT-kT /it7^t] , (41) 

sin(05i - ^2) [Pt-^t /it5it] ■ (42) 

The function /-j^ corresponds to the so-called Sivers effect Q . Here one usually assumes that such a function 
might arise due to initial state interactions and its contributions indeed do not depend on the lepton pair 
orientation as can be seen from the above two angular dependences. 



10 



Therefore, these structures are distinguishable from the T-odd asymmetry A^rp. However, it is important 
to note that Aj^j, can also effectively arise due to SM CP violation, hence this contribution must first be 
calculated before any conclusion about physics beyond the SM can be made. Also, this specific asymmetry will 
be accompanied by the product h'^{xi)hi{x2), thus it will suffer from the same drawback as Eq. (|^), namely that 
the transversity function for the antiquarks is presumably smaller than for the quarks, making this asymmetry 
hard (if not impossible) to detect at RHIC. But it illustrates how symmetry violation can be used in principle 
to disentangle SM asymmetries from new physics asymmetries. 



IV. CONCLUSIONS 



We have investigated a helicity non-flip double transverse spin asymmetry in vector boson production in 
hadron-hadron scattering, which was first considered by Ralston and Soper. It does not involve transversity 
functions and in principle also arises in W^-boson production for which we have presented the expressions. The 
asymmetry requires observing the transverse momentum of the vector boson, but it is not suppressed by explicit 
inverse powers of the large energy scale Q. However, as we have shown, inclusion of Sudakov factors suppresses 
the asymmetry at least by an order of magnitude compared to the tree level result. This suppression increases 
with energy approximately as a fractional power, numerically found to be a ~ 0.6. Moreover, the asymmetry 
is shown to be approximately proportional to xigi{xi) X2gi{x2), which gives rise to additional suppression at 
small values of the light cone momentum fractions. This implies that the asymmetry is negligible for Z and W 
production at RHIC and is mainly of interest at low energies (for 7* production). The strong suppression with 
respect to the tree level result will also have consequences for similar types of transverse momentum dependent 
azimuthal asymmetries in for instance e+e" annihilation at the Z mass scale, where the same strong suppression 
due to Sudakov factors is expected. 

We have also noted that unlike the transversity and CP-violating double transverse spin asymmetries, the 
helicity non-flip asymmetry Att(Qt) does not depend on the orientation of the transverse spin vectors compared 
to the lepton pair production plane orientation. This feature can be exploited to separate the different types of 
asymmetries. 



ACKNOWLEDGMENTS 



I would like to thank Les Bland, Gerry Bunce, Bob Jaffe, Rainer Jakob, Piet Mulders, Naohito Saito and 
Werner Vogelsang for helpful comments and discussions. Furthermore, I thank RIKEN, Brookhaven National 
Laboratory and the U.S. Department of Energy (contract DE-AC02-98CH10886) for providing the facilities 
essential for the completion of this work. 



[1] J. P. Ralston and D.E. Soper, Nucl. Phys. B152, 109 (1979). 

[2] F. Baldracchini, N.S. Craigie, V. Roberto and M. Socolovsky. Fortsch. Phys. 30, 505 (1981). 

[3] V.A. Korotkov and W.D. Nowak, [hep-ph/990849C| . 

[4] O. Martin and A. Schafer, Z. Phys. A 358, 429 (1997); 

V. Barone, T. Calarco and A. Drago, Phys. Rev. D 56, 527 (1997); 

O. Martin, A. Schafer, M. Stratmann and W. Vogelsang, ibid. 60, 117502 (1999). 
[5] R.D. Tangerman and P.J. Mulders, Phys. Rev. D 51, 3357 (1995); [hep-ph/9408305 . 
[6] C. Bourrely and J. Soffer, Nucl. Phys. B423, 329 (1994). 
[7] J.C. Collins, Nucl. Phys. B396, 161 (1993). 

[8] R.L. Jaffe and X. Ji, Phys. Rev. Lett. 67, 552 (1991); Nucl. Phys. B375, 527 (1992). 

[9] J.C. Collins and D.E. Soper, Phys. Rev. D 16, 2219 (1977). 
[10] D. Boer, P.J. Mulders and O.V. Teryaev, Phys. Rev. D 57, 3057 (1998). 
[11] D. Boer, Phys. Rev. D 60, 014012 (1999). 
[12] W. Vogelsang and A. Weber, Phys. Rev. D 48, 2073 (1993). 
[13] D. Boer, R. Jakob and P.J. Mulders, Nucl. Phys. B564, 471 (2000). 
[14] X. Ji, Phys. Lett. B 289, 137 (1992). 



11 



E. Mirkes and S. Willfahrt, Phys. Lett. B 414, 205 (1997). 

J.C. Collins, D.E. Soper and G. Sterman, Nucl. Phys. B250, 199 (1985); G.T. Bodwin, Phys. Rev. D 31, 2616 
(1985); 34, 3932(E) (1986). 

M. Ahmed and T. Gehrmann, Phys. Lett. B 46 5, 297 (1999). 

G. Bunce, N. Saito, J. Soffer and W. Vogelsang, [hep-ph/0007218 , 

R.L. Jaffe and N. Saito, Phys. Lett. B 382, 165 (1996). 

R. Jakob P.J. Mulders an d J. Rodrigues, Nucl. Phys. A626, 937 (1997). 

D. Boer, [iep-ph/991231l| . 

Abe et ai, Phys. Rev. Lett. 76, 587 (1996); Phys. Rev. D 58, 112003 (1998); 
Adams et ai, ibid. 56, 5330 (1997). 



E143 Collaboration, 
SMC Collaboration 



K. 
D. 



C. Balazs, Ph.D. Thesis, Michigan State University, hep-ph/9906422 



J.C. Collins, in Perturbative Quantum Chromodynamics, edited by A.H. Mueller (World Scientific, Singapore, 1989), 
p. 573. 

J.C. Collins, D.E. Soper and G. Sterman, Phys. Lett. 109B, 388 (1982); Nucl. Phys. B223, 381 (1983); 
Phys. Lett. 126B, 275 (1983). 

D. Boer and P.J. Mulders, Nucl. Phys. B569, 505 (2000). 
Davies and W.J. Stirling, Nucl. Phys. B244, 337 (1984). 
Weber, Nucl. Phys. B382, 63 (1992); B403, 545 (1993). 
Frixione, P. Nason and G. Ridolfi, Nucl. Phys. B542, 311 (1999). 
Boer, R. Jakob and P.J. Mulders, Nucl. Phys. B504, 345 (1997). 
Boer, R. Jakob and P.J. Mulders, Phys. Lett. B 424, 143 (1998). 



G.A. Ladinsky and C.-P. Yuan, Phys. Rev. D 50, R4239 (1994). 

A.V. Efremo v, O.G. Smirnov a and L.G. Tkachev, Nucl. Phys. B (Proc. Suppl.) 79, 554 (1999). 



V.L. Rykov, [hep-ex/990805C . 

D. Sivers, Phys. Rev. D 41, 83 (1990); 43, 261 (1991). 



12 



